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U.S.  Air  Force  Academy  (Department  of  Electrical  Engineering),  Air  Force 
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ANTENNA  PATTERN  SYNTHESIS  COMPUTER  PROGRAM 


1.  INTRODUCTION 

Optimization  methods  have  been  successfully  applied  to  the  pattern  syn- 

1 2 

thesis  problem  of  antenna  arrays.  ’ This  application  is  independent  of  the 
functional  dependence  of  the  array  parameters  and  therefore  treats  them  on  an 
equal  basis.  This  independence  makes  the  method  very  general  and  applicable 
to  a truly  large  class  of  problems.  Furthermore,  the  method  allows  the  intro- 
duction of  a variety  of  linear  and  nonlinear  constraints  in  the  synthesis, 
giving  the  designer  the  opportunity  to  simplify  the  antenna  construction.  The 

optimization  method  in  this  report  is  based  on  one  of  the  well  known  gradient 
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methods  called  PARTAN  ’ implemented  in  the  form  of  a computer  program.  The 
PARTAN  method  has  been  selected  due  to  its  simplicity  in  choosing  a step  size 
at  each  iteration  and  implementing  constraints. 

2.  THE  PARTAN  OPTIMIZATION  METHOD 

The  radiation  patter  for  0 = 0q  of  a planar  array  on  the  x-y  plane  can 
be  written  as 

M jk(x  cos<j)+y  sin<t>)sin0 

E(4>)  = l (AR  + j AI  ) e m m o (i) 

m=l  m m 

where  the  pair  (AR  , Al  ) specifies  the  real  and  imaginary  parts  of  the  mth 
m m 

current  element  and  the  pair  (x  , y ) determines  the  position  of  the  mth  ele- 

m m 

ment  in  the  x,y  plane.  M is  the  number  of  antenna  elements  and  k is  the  wave 
number  (k  = 2tt/A).  In  this  first  treatment  the  mutual  coupling  between 
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SUM C.WT  num> 


elements  is  ignored,  and  each  antenna  has  an  omnidirectional  pattern.  Later 
the  mutual  coupling  effect  is  considered  and  the  current  on  each  antenna  is 
dependent  on  all  other  antennas  relative  positions. 

Let  be  the  specified  or  desired  radiation  pattern.  The  error  over 

the  entire  synthesis  range  can  be  defined  as 


C1-  I »<VIE<*„>  - *„<♦»> 

n=l 


(2) 


e2  ' l "<V 

n=l 


lEC^I  - |E8(*n) 


(3) 


where  N is  the  number  of  <J>  directions  on  which  the  pattern  is  specified  and 
W($  ) is  a weighting  function  that  allows  the  synthesis  precision  to  be 
changed  over  certain  ranges  of  q>.  Here  W(<J>n>  is  a non-negative  number.  When 
is  used  for  the  error,  only  the  amplitude  of  the  specified  radiation  pat- 
tern is  required.  While  the  amplitude  and  phase  of  the  specified  radiation 
pattern  is  required  for  e . There  are  some  other  error  definitions,  but  only 
and  are  used  in  this  report. 

Now  the  problem  is  to  find  the  minimum  of  the  error  function  £ by  using 

the  PARTAN  method.  We  start  with  an  initial  guess  for  the  variables  (AR^,  AI^, 

m m 

Xm’  ^m^  ’ comPute  the  gradient  Ve^*  at  this  point  and  then  compute  the  new  values 

of  the  variables  (AR^,  AI^,  x^*,  y^)  as 

m m m m 


AR1  = AR°  - t |§r- 
m m dAR 


AI1  = AI°  - t 11=- 
m m dAI 


1 0 9e 

x = x - t ■» — 
mm  dx 


1 0 9e 

y = y - t -5— 
m m dy 


(4) 
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where  the  derivatives  are  evaluated  at  the  initial  point.  We  then  continue 
in  a similar  fashion  until  the  minimum  of  e is  found.  The  problem  is  to  find 
t,  the  step  size,  at  each  iteration. 

Using  a Taylor  series,  we  get  from  Eqs.  (1)  and  (4) 

. M n ~ 0 „ a 0 jk(x^cos<t>+y^sin0)sin0r. 

E <*)  m l { <ARm  " t a ar“>  + J <AI™  " t e 


m=l 


m 


3AI 


m 


3 0 a 0 

x {l  - jk  sin  0Q  cos  0 + -p-  sin  4>)t  + . . .}  (5) 

m ym 


This  is  used  to  calculate  the  error  function  e or  C ^ in  Eqs.  (2)  or  (3). 

After  Eq.  (5)  is  substituted  into  Eq.  (2)  or  (3),  and  the  higher  order  terms 
2 

than  t are  neglected,  we  get  a quadratic  equation  in  t for  or  Neglect- 


ing higher  order  terms  than  t is  equivalent  to  assuming  that 

a 0 aP° 

jk  sin  0^  (p—  cos  + g — sin  <J>)t|  « 1 
m ym 


(6) 


Thus  the  step  size  t of  each  iteration  is  found  in  close  form  from  the  quadra- 
tic equation  with  respect  to  t for  the  minimum  or  in  the  direction  Ve^. 
This  is  one  of  the  advantages  of  approximating  Eqs.  (2)  and  (3)  with  the 
quadratic  equation  of  t.  Even  if  Eq.  (6)  is  not  satisfied  in  the  initial 
iterations,  eventually  it  will  be  satisfied  and  the  minimum  or  can  be 
obtained . 


3.  APPLICATION  TO  SOME  SPECIFIC  SYNTHESIS  PROBLEMS 
3.1.  All  Antennas  Are  Excited  in  Linear  Antenna  Array. 

Let  us  now  apply  the  above  ideas  to  some  specific  cases. 
Let  us  calculate  by  using  Eq.  (5). 
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n=l 


1 N M 0 9e?  n 3e? 

i'J,  J,  (<AR.- t5Ar> + J(AI^- 3Ar>) 

m-i  m ra 

0 9ei  0 9ei 

jk{(xm  ■t^r)cosV(vt8r)sin4,n}sin0o 

■e  m E (<D  ) 

s n 


Neglecting  the  terms  of  order  higher  than  t , we  get 


l N 


n=l 

N 


a + b t + c t 
n n n 


= I ( (a  a ) + (a  b + ab)t+(ac  +ac  +b  b*)t^} 
n_^  nn  nn  nn  nn  nn  ntr 


where 


(7) 


(8) 


a = 
n " ' m 
m=l 


w 


t1  n n jk(x  cos<J>  +y  s.in<J)  )sin0A 

y (AR  + j AI°)e  m n m m 0 _ 

=1  m m 

u V rd£l  9£1  0 0 oti  dei 

bn  = {3aT  + 9 JKT  + -ik  Sin  0O  (ARra  + (cos  *n  9lT  + sil^n  W^)] 


jk(x°cos<f>  +y°sin<t>  )sin0A 
m Tn  m n 0 


M SE“  3£J  Se“  9e" 

cn  " „Vlk  1 Sar  + 1 1aF)<cos  ♦n  ST  + si“  jTr>  sl»  0 

m-  i m m m J n 


k2  sln  eo  ,.„0  . ...0  3E 


0 


ae?  2 


2 (ARm  + iAIm)[cos  4>n(^)  + sin  <Dn(^>]  } 


m 

^jk(x0cos$n+y0sin$n)sin0o 

Es  ^ n^  ’ 


K9) 


The  derivatives  are 
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3e*?  N , -jk(x°cos4>  +y°sin4>  )sin0 

if-  2-  I ReUE1^)  - Es(V)e  » » ” " 0 


3AR 


n=l 


9e 


3AI 


T-  = 2*  I iJtE1^)  - E8(<Dn))e 


-ik(x°cosd>  +y°sin<J)  )sin6f. 
J m n m n 0 


K!0) 


m 


n=l 


9e?  N . n n ~jk(x°cos4>  +y°sin<J>  )sin0 

^=2-  £ k cos  4)  I [ (E1(*  ) -E  (<J>  ))(AR°- JAI°)e  ] 

dx  nm  nsnmm 

n=l 


!!i 

9y„ 


0 N , ^ n -jk(x^cos<{)  +y^sin4>  )sin0„ 

l=2‘l  k sln  Im[(E  (0n) 'Es(4,n))(ARm"jAIm)e  " ^ 

n=l 


where  Re  is  the  real  part  operator  and  I is  the  imaginary  part  operator. 

Equation  (8)  was  applied  to  synthesize  the  specified  pattern  shown  in 
Figure  1,  which  has  a triangle  shape  and  has  a peak  at  <t  = 90°.  The  array 
elements  are  allowed  to  move  only  along  the  x-axis.  Therefore,  24  variables 
are  used  for  this  problem.  The  weighting  function  W was  set  to  1.  In  the 
initial  guess  all  amplitudes  have  been  set  equal  to  1 and  the  antennas  are 
equally  spaced  on  the  x-axis  a half  wavelength  apart. 

The  second  example  is  the  synthesis  of  the  pattern  of  Figure  2 with 
20  dB  side  lobe  level.  The  array  elements  are  allowed  to  move  only  along 
the  x-axis.  Therefore,  24  variables  are  used  for  this  problem  also.  The 
weighting  function  W provides  a convenient  tool  for  increasing  the  accuracy 
of  the  synthesis  in  certain  portions  of  the  desired  pattern.  In  this  example 
W is  equal  to  5 for  points  in  the  main  beam.  In  this  way  emphasis  has  been 
placed  on  obtaining  a pattern  which  has  the  correct  main  beam-width.  For  the 
side-lobe  region,  W was  set  to  either  0 or  1,  depending  on  whether 
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E (tp  ) - E(<p  ) is  positive  or  negative.  In  this  way,  points  below  the  speci- 
s n n 

fied  side  lobe  level  have  no  influence  in  the  computation  of  the  error 
In  the  initial  guess  all  amplitudes  have  been  set  equal  to  1 and  the  anten- 
nas are  equally  spaced  on  the  x-axis  a half  wavelength  apart.  The  final 
solution  is  shown  in  Fig.  2. 


3.2  Linear  Array  with  Only  One  Excited  Element  and  the  Others  Parasite 

The  currents  on  the  parasite  elements  are  a function  of  their  relative 
position  and  of  their  distance  to  the  excited  antennas.  The  initial  currents 
on  all  antennas  are  obtained  by  using  the  method  of  moments^,  and  from  this, 
equivalent  point  sources  are  obtained  when  the  constant  9-plane  radiation 
pattern  is  considered.  These  point  sources  correspond  to  AR^  and  jAI^  in 
Eq.  (1).  For  this  problem  let  in  Eq.  (3)  be  the  error  function,  since 
in  most  cases  we  are  not  interested  in  the  phase  of  the  radiation  pattern. 


4 ■ £ "(Vl 

n=l 


|E1(4>n) 


E (4>  ) 

s n 


(id 


The  derivatives  to  be  used  are 


i 

9x 


B 3|E°M>  )|  . 

■ l <Ie°<V 

n=l  m 


3e°  N 3 1 E°(<(>  )| 

3T  -2'  £-3TL-(|E  (V 

m n=l  m 


lEs(<f>n)l)W(^n) 


- |E  (<f  )|)W(4>  ) 
' s n n 


(12) 


The  radiation  pattern  is  from  Eq.  (1) 

M jk(x  cos<t>  +y  sin <p  )sin0 


E(<J>  ) = J (AR  + JAI  )e 
n L,  m m 

m=l 


m n m n 


0 


(13) 


where  AR  + jAI  is  the  function  of  the  positions  of  all  antennas, 
m J m 
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If  mutual  coupling  effects  of  all  antennas  are  considered,  the  equi- 
valent point  sources  have  to  be  obtained  after  each  iteration  by  using  the 
method  of  moments.  This  involves  a large  number  of  calculations  and  compu- 
ter execution  time.  One  way  to  avoid  this  problem  is  to  assume  that  only 
the  mutual  coupling  between  each  parasite  and  the  excited  antenna  is  impor- 
tant. Further,  we  assume  that  in  Eq.  (13) 

IC  -jk(r  -rC) 

AR  + JAI  = — e m m (14) 

m J m r 

(~f) 

r 

m 

c c 

where  I is  the  initial  equivalent  point  source  and  r is  the  initial  dis- 
ra  m 

tance  between  the  excited  antenna  and  the  mth  parasite.  For  the  excited 

antenna  it  is  assumed  that  (r  /rC)  = 1.  For  simplicity  let  us  assume  that 

mm 

the  excited  antenna  is  m = 1.  Then  we  get 


(x  -x  )2  + (y  -y.)2 
mi  mi 


v-  j<v*;>2  + <v^>2 

c c c c 

where  (x  , y ) and  (x, , y.. ) are  initial  positions  of  the  parasite  and  the 
m m l l 

excited  antenna  respectively. _ 

9EU(<}>n)  3E°(<fr  ) 

Now  the  derivatives  — t- and  ~ can  be  obtained  as 

dx  9y 

m 'm 


0 c c 

9E  (cj)  ) r -jk(r  -r  ) jk(x  cos$  +y  sin$  )-sin9„  , (x  -x, ) 

n Tc , ms  J m m J ' m Tn  7m  n' ° , ,..  . 1,  m 1 , 

— 5 = I ( — ) e e [jk  cos  - (jk+ — ) ] 

9x  m r n J r r 

n m mm 


for  m ^ 1 


9E°(4>  ) jk(x  cos(J)  +y  sin<£ 

n ,t  tc  i ni  nm  n 

— ~ = jk  I,  cos  e 

dx^  1 n 


? 


9E^(<J>  ) rC  -jk(r  -rC)  jk(x  cos<f>  +y  sin<)>  )-Sin.90 


n Tc , m x 

= ! ( )e 

9y  m r 

m m 


mm'-'  ' m 'n  'm  Tn'ri,  x , l^ym  yi^, 

e [jk  sin  ^ - (jk  +— ) — ] 

m m 
for  m ^ 1 


9E^(<}>  ) jk(x  cos<f  +y  sin<}>  ) • S i R Q„ 

n • i j i m rim  n 

— r = ikl,  sin  0 e 

9y.  1 Yn 


(16) 


Since  | E° (4>  ) J 2 = E°(<t>  )*E°(i  ) , we  get 


|E  (4>, 


9[EU($n)| 


n n 

3E°(4»n) 


3x 


■ Re  ‘-sr2-  E<V  1 


m 


Then 


9e°(4>_) 


9|E°(4>  )|  Re  [ 9x  E(<|,n)  ] 


3x 


ra 


|E°(4>n)| 


Similarly 

3E°(»  ) » 

3|E°(V|  Re  <-^T-  E<V  1 


9y„ 


E%n)| 


> (17) 


Substituting  Eqs.  (16)  and  (17)  into  Eqs.  (12),  we  can  obtain  the  value  of 
the  derivatives  to  start  the  iteration. 


From  Eqs.  (4),  (13),  (14)  and  (15),  we  get 


EV  - I'e 


M 


jk[  (x^-t  (902/9x1))cos$n+(y^-t  Oe^/Sy^sin^JsinGQ 


r Tc  c me 
+ ) Ire 

L~  mm 
m=2 


jkrC  -jk  >/(: Vc  (3e2/3xm)“xl) 2+(: Vc  (8e2/3ym)~yl)  2 


& 


tOe^/dx  ) - *°)2  + (y°  - tOe^/Sy  ) - y?)2 
z m l m Z m 1 


jk[  (x°-t  (9e!?/9x  ))cos0  + (y^-t(9e^/9y  ))sin<t>  ]sin0 
m z m n m Z m no 


(18) 
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By  using  a Taylor  series  and  neglecting  the  higher  order  terms  than  t‘ 
get 

eV)  = l V1  + Blt  + B2t2) 

m=l 

where 


c o c 0 0 

r -ik(r  -r  ) jk(x  cos  <p  +y  sin<{>  )sin0„ 
. Tc,  m , J m m m Tn  m n 0 

A = I ( ) e e 

m m o 
r 

m 


B,  = 


- jk(g  cos  <b  + g sin  d>  )sin6n  +D, 
J x n y n 0 1 

f 0 


D,  = 


(xm  ~ Xl)8x  + (ym  ~ yl)g 


m 


1 (2k  + -V 

r 


m 


m = 1 


m £ 1 


B = - 


2 2 

k sin  0Q  2 

2 (gx  COS  ^n  + 8y  sin  V + °2 


D„  = 


0 0, 


0 0,  2 


o t 


t 


_S_ 
(rm)2 


m = 1 


2 2 

+ g. 


f<V’‘l>g*  + <Vy-^}  (3  + 3k r»)  - a + 3kr»)] 


(tV 

m 


, 0 0.  .,  0 0, 

(x  -x,  )g  + (y  -y,)g  / 

m l/sx  VJm  ,1  ey 

Jk  — — • < g cos  <p  -*■  g s]  n (j> 

a n y n 


, 0.2 

(O 

m 


, 0 0,  ,,  0 0, 

(VX1)8X  + (Vyl)8v  ] 

f n * 1 

r J 


3e; 


8„  = 


x 9x 


9e 


8.,  = 


y 9y„ 


, we 
(19) 
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Eq.  (19)  can  be  expressed  as 
E^(<J>  ) = a + bt  + ct2 


' n 
where 


M 

a = l A , 
m-1  m 


M M 

b = y A B,,  and  0=7*  A B„ 
, m 1 m ^ 


m=l 


m-1 


(21) 


Then  Eq.  (21)  is  substituted  into  Eq.  (11)  and  we  have 
N 


hi  " l < J< 


n=l 


^ a o o 

a + bt  + ct  )(a  + b t + c t ) - lEs(<J>  )|) 


(22) 


Neglecting  the  terms  of  orders  higher  than  t“ , we  have 

1 r . I * * * * * 2 , ‘ 

e?  = l ( Jaa  + (ab  + a b)t  + (ac  + a c + bb  )t  - |E  (<J>  )|) 

n=l  s n 


I [ fiT-  |Ea(»  )|  + t + (^"  +7a  c, + bb"  - (a-b>^;^-)t2] 

n=l  ‘ 2 J aa  2J  aa  4( J aa  )J 


* * 


2 2 


- I I E (0  ) I ) 2 + l l ab_^a  b ( JI 7 - |E  (<{»  ) 1 ) 1 1 

n=l  S n n=l  Jaa*  S n 


N 

l 

n=l 


s n 


+ l I ((  I E («J|>(— -jg=v'b->  + 1*  — | E <4  >|}]t‘  (23) 


aa 


* * ) 

> 

.3  ' s ' ' n' 


4(  Jaa*)' 


Equation  (23)  is  a quadratic  function  of  t and  the  step  size  t is  found 
easily  for  the  minimum  for  This  process  continues  until  the  change  of 

is  small  enough  at  each  iteration. 

The  assumption  of  Eq.  (14)  is  now  under  verification.  Preliminary 
results  have  shown  it  to  be  solid  if  in  computing  the  current  in  the  parame- 
ters the  presence  of  the  two  adjacent  antennas  are  considered. 
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4. 


CONCLUSIONS 


The  computer  programs  for  the  cases  discussed  in  Section  3 have  been 
written  and  debugged.  Once  the  validity  of  Eq.  (14)  or  any  equivalent  formu- 
lation is  verified,  we  can  proceed  with  the  method. 
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FIGURE 


METRIC  SYSTEM 


BASE  UNITS: 


Quantity 

_ Un|t 

SI  Symbol 

Formula 

length 

metre 

m 

mass 

kilogram 

kg 

time 

second 

s 

electric  current 

ampere 

A 

thermodynamic  temperature 

kelvin 

K 

amount  of  substance 

mole 

mol 

luminous  intensity 

candela 

cd 

SUPPLEMENTARY  UNITS; 

plane  angle 

radian 

rad 

solid  angle 

steradian 

tr 

DERIVED  UNITS; 

Acceleration 

metre  per  second  squared 

rrv* 

activity  (of  a radioactive  source) 

disintegration  per  second 

(disintegration  1/s 

angular  acceleration 

radian  per  second  squared 

rad/s 

angular  velocity 

radian  per  second 

rad/s 

area 

square  metre 

m 

density 

kilogram  per  cubic  metre 

kgm 

electric  capacitance 

farad 

F 

A-s/V 

electrical  conductance 

siemens 

S 

A/V 

electric  field  strength 

volt  per  metre 

Vm 

electric  inductance 

henry 

H 

V-a/A 

electric  potential  difference 

volt 

V 

WA 

electric  resistance 

ohm 

V/A 

electromotive  force 

volt 

V 

W/A 

energy 

joule 

1 

N-m 

entropy 

joule  per  kelvin 

|/k 

force 

newton 

N 

kg-m/s 

frequency 

hertr 

Hi 

(cycle  Vs 

illuminance 

lux 

lx 

1mm 

luminance 

candela  per  square  metre 

cd/m 

luminous  flux 

lumen 

im 

cd-ar 

magnetic  field  strength 

ampere  per  metre 

A/m 

magnetic  flux 

weber 

Wb 

V-s 

magnetic  flux  density 

tesla 

T 

Wbim 

magnetomotive  force 

ampere 

A 

power 

watt 

W 

)'* 

pressure 

pascal 

Pa 

N/m 

quantity  of  electricity 

coulomb 

C 

A-s 

quantity  of  heat 

joule 

1 

N-m 

radiant  intensity 

watt  per  steradian 

Wsr 

specific  heat 

loule  per  kilogram-kelvin 

I'kg-h 

stress 

pascal 

Pa 

N/m 

thermal  conductivity 

watt  per  metre  kelvin 

Wm-h 

velocity 

metre  per  second 

ms 

viscosity,  dynamic 

pascal-second 

Pa-s 

viscosity,  kinematic 

square  metre  per  second 

m/s 

voltage 

volt 

V 

W/A 

volume 

cubic  metre 

m 

wavenumber 

reciprocal  metre 

(wavel/m 

work 

joule 

1 

N-m 

SI  PREFIXES 

Multiplication  factors 

Prefix 

SI  Symbol 

1 000  000  000  000  !()'> 

Inra 

T 

1 000  1X10  000  = to* 

gig* 

<: 

1 000  000  --  10* 

mega 

M 

1 000  - 10' 

kilo 

k 

100  * 10' 

hecto* 

h 

10  - 101 

dekg* 

da 

0 1 = 10-’ 

decl" 

d 

001  - io-» 

1 entl* 

c 

0 001  * 10*» 

mill! 

m 

ii  ooo  ooi  itr* 

micro 

M 

0 000  000  001  10- • 

nano 

n 

o ooo  ooo  non  ooi  • 10  '» 

pirn 

n 

0 000  000  (MM)  0OO  001  10  '« 

inmto 

F 

0 000  000  000  000  000  001  10  '• 

alio 

M 

* To  be  avoided  where  poHible 

•US  GOVONMfWI  MlNtlXf,  Orrici  1*7-  T|«  OJJ  |«t 


MISSION 


Rome  Air  Development  Center 


RADC  plans  and  conducts  research,  exploratory  and  advanced 
development  programs  in  command,  control,  and  coammnicatlot 
(C3)  activities , and  in  the  C3  areas  of  informatior,  science 
and  intelligence.  The  principal  technical  mirslon  areas 
are  communications,  electromagnetic  guidance  ind  control, 
surveillance  of  ground  and  aerospace  objects,  intelligence 
data  collection  and  handling,  information  system  technology 
ionospheric  propagation,  solid  state  sciences,  micromave 
physics  and  electronic  reliability,  maintainability  and 
compatibility . 


